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Abstract. Relativistic theories of nuclear matter are discussed in a new pespective. First the chiral 
character of the scalar nuclear field is introduced in the framework of the linear sigma model. With 
the assumption that the nucleon mass originates in part from the coupling to the quark condensate 
it is possible to relate the optical potential for the propagation of the scalar field to the QCD scalar 
susceptibility of the nucleon, on which indications exist from the lattice evolution of the nucleon 
mass with the quark mass. Constraining the parameters of the nuclear scalar potential by the lattice 
expansion parameters a successful description of the nuclear saturation properties can be reached. 

^ \ Keywords: Chiral symmetry. Confinement, Nuclear Matter Saturation 

In this meeting devoted to the scalar mesons it is appropriate to talk of nuclear physics 
^ ■ and I am one speaker in charge of this task. Indeed it is clear that the existence of a scalar 

rn \ meson coupled to nucleons has big implications for nuclear physics. This idea has been 

exploited in the relativistic models of the nuclear binding of Walecka and Serot [111]. 
Here the binding results from the balance between the attractive sigma exchange and 

t;;j^ . the repulsive omega exchange. This model had impressive successes, in particular for 

O \ what concerns the spin-orbit coupling. The main new idea which has emerged since 

the introduction of this model is the concept of a nucleonic response to the scalar field 
proposed by Guichon [0,|3|]. The nucleon is not inert but it responds to the presence of 
the field. This concept in itself is not new and occurs for instance in pion propagation. 

/\ \ The original aspect of the quark meson coupling model (QMC) is that it implies the 

j^ ■ quark struture of the nucleon. We will come back to this model. Here we propose the 

existence of a link between the scalar parameters of the scalar potential and the QCD 
lattice data iJ] and I will try to convince you that it is natural and to be expected. 

If one talks about a link between nuclear physics and QCD, it is natural to invoke 
effective theories which simulate QCD at low energies in terms of hadronic degrees of 
freedom and we use as a starting point the linear sigma model. Here chiral symmetry 
is implemented with the introduction of two fields, the o, scalar isoscalar, and the K, 
pseudoscalar isovector, which are chiral partners. In this model the part of the lagrangian 
which breaks chiral symmetry, which in QCD is -^smbreak ~ ~2m^(7^, is replaced by 
■^vvmbreak ~ f^'^K ^- ^^^ ^^ ^^^ spontancous Symmetry breaking, the expectation values 
of the order parameter are non vanisishing, with the vacuum values {qq)vac in QCD 
and (a) = f^- The equivalence of the symmetry breaking parts of the QCD and model 
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lagrangians implies the following equivalence between the quark density and the scalar 

field: 

qq{x) o{x) 

{qq)vac h 

Taking expectation values we find that the quark condensate in the nuclear medium is 
governed by the expectation value of the sigma field. 

We have now to incorporate the concept of the nuclear scalar field in this theory. A 
natural and simple identification would be that of the nuclear mean scalar field with the 
expectation value 5 {o) = {o) — fn-lt would make life simple because a unique quantity 
8{o) would govern the evolutions of the condensate and of the effective nucleon mass 
in the nuclear medium, A/^, with : 

M* {qqjp)) ^ ^ 5{a) ^ ^ Syyp;^ ^^^ 

Mn {qq)vac fn f^ml 

where p^ is the nucleon scalar density and Ea^ the nucleon sigma commutator. As 
the condensate evolution is known for independant nucleons from the value of Eyy, 
as written above, the mean scalar field would be known to first order in density, in a 
model independent way, with a value of ^ 30MeV at normal density. Unfortunately 
this simple identification is not legitimate, as emphasized by Birse [5]. The pionic 
contribution to the nucleon sigma commutator contains non analytical terms in m^. It 
follows that the NN interaction which governs the mass evolution would have terms 
of order m;^, which is not allowed by the chiral constraints. Of course it is always 
possible to add other exchanges, involving in particular two-pion states, in such a way 
that the unwanted terms cancel. But this is a very cumbersome procedure which is not 
practiced. We have proposed |S] a short cut with the suggestion that it is interesting to 
go from the cartesian coordinates, o and n, to polar coordinates, S and according to : 

o + ix-Tt = SU = {f^ + s) exp ( it ■ 4/1% ) . where S is associated with the radius of the 

chiral circle and (^ / f^ the polar angle. By doing so we go from a linear representation 
to a non linear one. Contrary to the usual procedure of the non linear model we do not 
freeze the radius of the chiral circle but we allow for a change in the nuclear medium, 
where it goes from its vacuum value, f^, to f-^ + s. We have suggested to link the 
nuclear field to this radial mode with the identification of the background mean nuclear 
scalar field with the mean value s. Since we are in a non-linear representation all chiral 
constraints are automatically satisfied. We point out that the bare masses of the o and the 
s are identical but these two fields differ by their off- shell couplings to two-pion states. 
In particular, due its derivative couplings, the s field decouples from low energy pions. 

Note that the extreme simplicity reached in the previous identification is lost. We 
have now two scalars, one 5(a) which is not invariant under chiral transformations 
governs the condensate evolution, and s which is a chiral invariant, governs the nucleon 
mass one. However they are not unrelated. Expressing the sigma field in terms of the 
polar coordinates we obtain for the condensate evolution, keeping only the lowest order 
terms : 

Mp))=i + jL_M. (2) 

{qq)vac fn 2/2 ■ 



Here the quantity ^0^) represents the scalar density of virtual pions in the nuclear cloud 
divided by the pion mass. For the nucleon mass instead : 

-^ = 1 + ir. (3) 

The nuclear scalar mean field still influences the condensate. However in order to extract 
it from the condensate, the effect of the pion cloud has to be separated out, which 
introduces an unavoidable model dependence. 

So we are now at peace with chiral symmetry but we face another problem. This chiral 
invariant nuclear scalar field is subject to the chiral dynamics. The mexican hat potential 
of the chiral theory contains a three scalar coupling term which produces a lowering of 
the sigma mass in the medium : 

^a ^m^- -tPn^ (4) 

In 

or equivalently for the existence of a tadpole crA'^ amplitude : 

ToN = -^8s/f7f (5) 

Here gs is the a-nucleon coupling constant. This lowering is a large effect, about 30% 
decrease of the mass at normal density. It means more attraction with increasing density, 
such that a collapse occurs, instead of saturation, as pointed out long ago by Kerman 
and Miller [7]. The o model simply does not provide a viable theory of nuclear matter. 
Phenomenologically, the problem can be cured with the introduction of a response of 
the nucleon to the scalar field, in the form of a polarizability coefficient, K^r such that the 
nucleon mass evolves, as in QMC, according to : 

1 , 

M^=Mn+ gss + - Kns\ (6) 

With a positive value of the response, repulsion is introduced and it is possible to reach 
saturation. For us [8], with an empirical value of Kpj such that approximately two thirds 
of the tadpole effect are cancelled, saturation is possible. But phenomenology is not our 
aim here and no link with QCD emerges. The link will appear with the study of the 
QCD scalar polarizability, Xs^ defined as Xs = d{qq) /dniq, the derivative of the order 
parameter {qq) with respect to the bare mass parameter m^ which breaks explicitely 
chiral symmetry. For the nuclear medium we define the susceptibility Xs ^^ ^'J'^h a way 
that the vacuum value is subtracted off, and only the condensate evolution enters. The 
susceptibility represents also the propagator of the fluctuations of the order parameter, 
d{qq){x), which in the sigma model is simulated by 5o{x). The susceptibility is thus 
expressed in the model in terms of the o propagator, with : 



Xs - ^ 72— 

J n 



1 1 



m^2 ^2 



(7) 



where m^ is the in-medium sigma mass, as modified by the tadpole term. Expanding 
the r.h.s. of eq. |7]to first order in the nucleon density we obtain : 

The meaning of this linear term in density is clear. We are dealing with the modifications 
brought into the vacuum susceptibility by the presence of the nucleons. Among these, 
part arises from the susceptibility of the nucleons which individually respond to a 
modification of the quark mass. The factor in front of the density therefore represents 
the scalar susceptibility of a free nucleon which is found as [|9|] : 

The linear sigma model thus predicts the existence of a negative component in Xs which 
is associated with the scalar field and linked to the tadpole amplitude. After having 
obtained this result it was possible to derive it directly in the model, as the derivative of 
the nucleon scalar charge with respect to the symmetry breaking parameter but in this 
case the connection to the tadpole amplitude does not emerge. 

This is the sigma model prediction. In real life, is there any evidence for this compo- 
nent ? Possibly, in the lattice simulations of the evolution of the nucleon mass with the 
quark mass (equivalently the squared pion mass) llO, [ill]- These results do not cover 
the physical mass but only a region of the quark mass above ^ 50MeV with several 
data points above this value. In order to extract the physical nucleon mass, an extrapola- 
tion has to be performed. For us the fact of having several data points at varying quark 
masses is interesting because the successive derivatives of the nucleon mass with respect 
to the quark mass provide the nucleon sigma commutator and the scalar susceptibility, 
according to : 

2 dM 
LN = m^^-^, (10) 

and : 

fn dml \mlJ 

These however are total values which include the pionic contributions, which should 
be removed for our nuclear physics purpose. Fortunately the pion loop contribution to 
the nucleon mass evolution has been separated out with great care in the works of ref. 
nid. Inn . The reason is that it contains non analytical terms in the quark mass which 
prevent a small mass expansion. The pion contribution depends on the nN form factor. 
Different form factors have been used with an adjustable cutoff parameter A, which is 
fitted. The non-pionic mass is expanded in niq (equivalently in m|) as follows : 

MN{ml.) =ao + azml +a4m^+ E;r(m;j,A). (12) 

The best fit values of the parameters a4 and a2 show little sensitivity to the shape of the 
form factor, witha4 ^ —O.SGeV^^, while a2 — l.SGeV^^ [ | lli] . From the expansion of 



eq. (fT2l) we deduce : 

Y^'^^'^-P'"'' = 2mqQs = ml—^=a2ml + 2a4ml ~a2m| = 29MeV. (13) 

It is largely dominated by the a2 term. The quantity Qs defined above is the total scalar 
quark number of the nucleon, that we denote also as the nucleon scalar charge. In turn 
the nucleon susceptibility is : 

^N,non-p,on ^^ (gg)L <^^ /g^n^\ ^ (^^^^ ^ SAGeV-' (14) 

The non-pionic susceptibility is found with a negative sign, as expected from the sigma 
model. If signs are right, are the magnitudes also compatible with the sigma model ? In 
this case the scalar charge is : 

^ {qq)vacgs 

As for the suceptibility the identification of the sigma model result and the lattice one 
leads to: 

^^Myon-,.on ^ 2(&£ /-3^\ ^ (^^^^ ^^^^ 

gs \ f^ J H 

which, using the relation (fT3l) between Qs and a2 gives : 

-^4-^ ^ V I c.l^ = 3.5GeV-\ (16) 

2 gsfnmn^ 2 gsfn 

seven times larger than the lattice value, —a^ = 0.5 GeV^-^. The sigma model is contra- 
dicted also by the value of the expansion parameter a4, a new failure for this model ? 
No, it is in fact the same as the previous one since the susceptibility and the tadpole 
amplitude are proportional and the sigma model prediction fails for both, as it should. 
This full consistency between the nuclear physics and the QCD lattice results is even 
gratifying. The question is if we can go beyond this failure. There is a need for a com- 
pensating terms. The common cure is found in confinement, which introduces a positive 
response of the nucleon to the scalar field, as in QMC. This model is a pure bag model in 
which the nucleon mass entirely arises from confinement. Moreover the chiral caracter 
of the nuclear scalar field is not discussed. Here instead we will assume that the nucleon 
mass originates in part from confinement and in part from the condensate. In addition 
we keep the assumption that the nuclear scalar field is the chiral invariant field discussed 
previously which influences the quark condensate. 

The first question is : in this mixed situation what is the relation between the nucleon 
scalar susceptibility and the oN scattering amplitude ? We will illustrate this point in 



a hybrid model of the nucleon, similar to the one introduced by Shen and Toki [11211 
which consists in the following. Three constituant quarks, described by the Nambu- 
Jona-Lasinio model, move in a non-perturbative vacuum, their mass M is generated 
from the coupling to the quark condensate. They are kept together by a central harmonic 



potential which mimics confinement and the effect of the color string tension. Although 
oversimplified the model provides an intuitive picture of the role played by confinement. 
The nucleon mass, because of the confining force, becomes 3E{M), the M dependence 
being fixed by the type of force. For illustration we take a harmonic potential of the form 
{K/4){1 + 7o) r^, which leads to the expression : 



Mn = 3E = 3\M+-J-——-]. (17) 




It is increased as compared to the value, 3M, for three independent quarks. With our 
assumptions the presence of the mean nuclear scalar field in the medium modifies the 
condensate and hence affects the mass M. The derivative dM/ds = gq — M/f^ has the 
non- vanishing value of the NIL model. The scalar coupling constant of the nucleon is : 



dM^ _ dE_ 



^5 = ^f = 3^,^. (18) 



The nucleon scalar charge, Qs, writes : 

3 dE 3 dE dM 
2 oniq 2 dm oniq 

where 3dE/dM = 3 {cs = 3E + 3M/3E+M) is the scalar number of constituant quarks. 
As E > M,cs < 1 , it is reduced as compared to a collection of three independent quarks. 
The nucleon scalar susceptibility, Xs ^ given by the next derivative, is composed of two 
terms arising respectively from the derivative of cs and from that of dM/dmq : 



^^ dm,, 2 



dcsfdMy d^M 



dcs 24(^2 _m2) 
with 3t| = ^ ^. (20) 



Notice that this last derivative vanishes in the absence of confining force, when E = M, 
and that it is positive since E > M. Therefore the first part of the expression of Xs 
represents the part of the susceptibility originating in confinement which is positive as 
in QMC. While instead the other term, proportional to the susceptibility of a constituant 
quark, d^M / d^rn^, is negative. Thus a compensating effect is possible. 

Is there a similar compensation in the crA^ scattering amplitude ? Two terms as 
well contribute to TaN- One is the tadpole term on the constituant quarks. For each 
constituant quark the tadpole amplitude is t(,i^ = —3gq/fjz. Multiplying by the scalar 
number of constituant quarks which is 3cs, the tadpole amplitude for the nucleon 
writes Ta^"^^"^^ = —3gs/ fit, the same expression as in the linear sigma model. We 
can compare this amplitude to the part of the nucleon susceptibility arising from the 
constituant quark susceptibility, second term on the r.h.s. of the expression (eq. [20l) of 
X^ . The constituant quarks as described by the NIL model obey at the quark level, the 
same relations as in the linear sigma mode. In particular the relation between the quark 
susceptibility and the tadpole amplitude is : 

2 



d^M _^(dMy f-3gq\ 



Multiplying both members by 3cs we find that the nucleonic suceptibility originating 
from that of the constituant quarks one is related to the tadpole amplitude by the same 
ratio 2 Qg/g^ as previously. As for the part originating in confinement, the amplitude K>/ 
is obtained as the second derivative of the nucleon mass with respect to the scalar field : 

d^E dcs fduV 

The ratio, r„,, between the part of the nucleon scalar susceptibility due to confinement 
and fCyv is 

the same ratio r as was previously found. Adding then the two components we obtain : 

N _ ^Qs rptotal /^ox 

AS — 2~ crA' ^ ^^ 

Ss 

The introduction of confinement has preserved the relation between the scalar suscep- 
tibility and the oN amplitude. Numerically our model is not successful. It does not 
produce enough cancellation of the tadpole effect but it is conceptually important to un- 
derstand the role of confinement. The relation between the susceptibility and the crA'^ 
amplitude established in the framework of this model is actually independent of the par- 
ticular form of the function ^(M) . A more sophisticated model of the nucleon is needed 
for a quantitative description. 

The relation (l23l) allows to go to the next step of our approach which consists in using 
the lattice data to fix or at least constrain the nuclear physics parameters of the scalar 
interaction part. In this approach for instance the medium effects in the propagation of 
the scalar field directly follow from the eq.(l23l) and we can write the o propagator as : 



-{Dl)-'=ml + gf-^p^==ml-0.5fpl (24) 

CI2 In 

The sigma mass is considerably stabilized with respect to the pure sigma model, where 
the coefficient in front of the last factor was the tadpole one, i.e., 3 instead of 0.5. 
The cancellation is so large that it leads to the question if we could not altogether 
forget the chiral effects as well as the nucleonic response. The answer is definitely 
no. The sigma propagator is not the whole story and there are important three-body 
forces which involve a different combination, with less cancellation, as shown below. 
In the expression Q for the nucleon mass evolution we perform a field transformation 
introducing a new scalar field, u = s + {kj^s^ /2gs), such that only a linear term in u 
enters this evolution. Expressed in term of the u field, the chiral mexican hat potential 
takes the form : 

y*'-''' = y = ^ {s^+j + ..\ = ^ ("' + j(i -2c) + ..)\ . (25) 



where we introduce the dimensionless parameter C = {KNfjt)/{2gs). In the formulation 
with the u field the three-body forces are totally contained in the u^ term of this potential : 

2 -3 
y three-body ^^^^i_2C). (26) 

2 fit 

As a < 0, this force is repulsive for C > 1/2, which is actually the case in our fit. 
Without confinement, i.e., for C = 0, the chiral potential alone leads to attractive three- 
body forces. The balance between the effects of the chiral potential and of the nucleonic 
response is not the same in the propagation of the scalar field and in the three-body 
forces. In the first case the amplitude T^*^"' = 'igs/fn + K"^ = {'igs/fn) (1 - 2C/3), 
while in the three-body forces the combination is 1 — 2 C. With C of the order one, 
which is the value in our fit, the cancellation of the tadpole term is nearly complete in 
jtjnai^ while there is an overcompensation in the three-body potential which becomes 
repulsive. The existence of repulsive three-body forces in relativistic theories is strongly 
supported by the nuclear phenomenology 11131 Il4l1 . They play an important role in the 
saturation. 

Going beyond the Hartree approximation where the pion does not contribute we have 
introduced [15] its contribution which occurs via the Fock term and the correlation 
one (including A excitations). For consistency we have also introduced the rho meson 
exchange. With the introduction of the pion (or rho) exchange interaction it is necessary 
to consider also the short range interaction in the spin-isospin channels in the form of 
contact interactions governed by the Landau Migdal parameters, ^/vAf,.?A?A'.?AA' where 
the indices refer to the type of hole or particle. 

Our fit which leads to a successfuU description of the nuclear binding goes as follows. 

a) For the quantity gs/fn^ we take the lattice value: gs/fn^ — 0.1/ fn = ISGeV^^. 
This corresponds, to leading order in density, to a mean scalar field of 20MeV at normal 
density. We also need separately the scalar coupling constant ^5, for which we keep the 
value of the linear sigma model, gs = Mf^/ f^^ , although confinement may introduce 
some deviation. The corresponding o mass is then ~ 800 MeV. 

b) The quantity C = {K]^fji)/{2gs) which determines the nucleonic response is al- 
lowed to vary near the lattice value, Ciattice = 1-25. The fit gives a value which is close, 
C=\. 

c) For the vector part of the potential, the omega mass has the experimental value, 
mco = 7S3MeV, while the coupling constant g^ is a free parameter. The fit value, gcj — 7, 
is close to the vector dominance quark model one, g^g ~ 8. 

d) We take for the nN form factor a dipole form with A = 0.9SGeV. It leads to a 
pionic sigma commutator a^""^ = 21MeV which, added to the non pionic value from 
the lattice yJ^"~P"'" — 29MeV, gives a total value E^v = 50MeV, in agreement with the 
experimental value. 

e) For the Landau-Migdal parameters of the spin-isospin interaction we use the latest 
information from spin-isospin physics [16] : ^^^ = 0.1 ,g'^^ = 0.3,^^^ = 0.5, producing 
a large suppression of the correlation terms. 

In summary we have found a full consistency between lattice QCD and the binding 
properties, relying on the following assumptions. The chiral invariant nuclear scalar 
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(a) (b) ('^) (d) 

FIGURE 1. Origin of the sigma mass in the vaccuum from ref. [17] (la and lb) and its modifications 
in the nuclear medium (Ic and Id). The quark lines (continuous lines) represent constituent quarks. 



field affects the quark condensate. The nucleon mass has a mixed origin, in part from 
confinement and in part from the quark condensate. With this limited set of assumptions 
it is possible to fix or constrain the scalar parameters by the lattice expansion ones, 
leading to a successful description of the nuclear binding properties. Confinement is an 
essential ingredient to reach saturation. It limits the attractive effect of the chiral mexican 
hat potential. The consistency between lattice data and the saturation properties, which 
is not a priori acquired, confirms the picture that the scalar potential which binds the 
nucleus reflects the modification of the QCD vacuum in the nuclear medium. This is true 
not only at the level of the mean scalar field related to the quark condensate evolution, 
pionic component removed, but also at the level of the three body forces, which are 
governed by the nucleonic QCD scalar susceptibility. 

This consistency between QCD and the nuclear potential is in itself remarkable but 
one may wonder wether there is a deeper necessity and significance for the introduction 
of both the tadpole term and the confinement one in the sigma propagation. After 
preparation of this talk I read an article by Delbourgo and Scadron [17] which may 
contain the answer to this question. In this work they introduce the linear sigma model 
at the quark level. Their lagrangian has a simple form which does not include the 
mass terms for pion and sigma, nor the mexican hat potential. These are generated 
dynamically by quark loops. In particular for the sigma the two graphs which build 
the o mass are the tadpole term involving a quark loop and the creation of a quark and 
antiquark pair by the o (fig. la and lb). These are constituent quarks whose masses obey 
a gap equation. This is the vacuum situation. Consider now how these two contributions 
are modified in the nuclear medium. The first one is changed by the extra tadpole term 
from the constituent quarks of the nucleons (fig. Ic). As for the second one, the creation 
of a quark antiquark pair by the sigma field, through the interaction of the sigma with 
a constituant quark of the nucleon, as shown in the graph Id, represents precisely the 
nucleonic response arising from confinement which comes from Z graphs. These two 
terms have opposite signs. This shows that the two contributions that we have considered 
in the sigma propagation in nuclei go together and are necessary as they represent 
the medium modifications of the two contributions which build the sigma mass in the 
vacuum. 
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